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ON  THE  STABILITY  OF  VORTEX  MOTIONS  IN 
COMPRESSIBLE  STRATIFIED  FLUIDS 

INTRODUCTION 

Vortex  shedding  behind  an  axisymmetric  body  submerged  in  stratified  fluids  has  been  puzzling 
many  researchers  because  of  the  intrinsic  vortex  pattern  developed  in  the  late  wake.  When  an  axisym¬ 
metric  body  is  towed  through  a  stratified  fluid  at  a  particular  Reynolds  number,  vortices  are  first  shed 
three  dimensionally.  However,  the  gravitational  effects  induced  by  the  density  stratification  soon  inhibit 
the  vertical  motion.  The  resultant  vortex  structure  is  vertically  oriented  and  resembles  the  two- 
dimensional  Kirman  vortex  street  behind  a  bluff  body  if  observed  from  the  gravitational  direction. 
Even  though  the  mechanism  behind  the  development  of  this  vortex  structure  is  still  not  clear,  the  inhi¬ 
bition  of  the  vertical  motion  can  qualitatively  be  viewed  by  examining  the  motion  of  fluid  particles  in 
the  gravitational  force  field.  When  a  sphere  or  other  axisymmetric  body  is  towed  through  a  stratified 
fluid,  the  fluid  particles  in  four  locations  are  of  particular  interest.  Figures  la  and  lb  respectively  show 
the  top  view  and  side  view  of  the  sphere  and  of  the  four  fluid  particles  A,  B,  C  and  D  being  considered. 
Since  the  fluid  is  stably  stratified  in  the  vertical  direction,  the  densities  of  particles  A  and  B  are  equal 
while  the  density  of  particle  C  is  lighter  than  that  of  particle  D.  To  create  a  rotational  motion  at  the 
onset  of  shedding  means  that  either  particles  A  and  B  or  particles  C  and  D  have  to  interchange  their 
positions.  The  interchange  of  particles  A  and  B  requires  no  work  done  in  the  gravitational  force  field. 
The  interchange  of  particles  C  and  D,  however,  requires  a  work  done  equal  to  the  increase  of  the 
potential  energy  at  the  new  location.  Accordingly,  nature  will  take  an  easy  way.  With  the  vertical 
motion  suppressed  by  the  gravitational  forces,  the  resultant  vortex  motion  will  be  confined  in  a  rela¬ 
tively  two-dimensional  motion  reminiscent  of  the  Karman  vortex  street  behind  an  infinitively  long 
cylinder  when  viewed  in  the  vertical  direction. 

Manuscript  submitted  on  June  18,  1981 
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Fig.  lb  —  Side  view  of  a  sphere  towed 
in  a  stably  stratified  fluid. 


9 


To  understand  the  behavior  of  such  vortex  motions,  we  need  to  mathematically  consider  a  more 
general  type  of  flow  profiles  which  can  be  used  to  describe  this  kind  of  flow  behavior.  Most  of  the 
steady  profiles  considered  in  parallel  flows  vary  only  in  one  direction,  i.e.,  they  are  single  coordinate- 
dependent,  even  though  they  may  possess  more  than  one  independent  velocity  component.  Relatively 
little  attention  has  been  given  to  flows  varying  in  more  than  one  direction.  The  vortex  motion  being 
considered  is  of  double  coordinate-dependency.  A  constraint  relation,  required  to  satisfy  the  pressure 
balance  condition  at  all  points  within  the  flow  domain,  may  play  an  important  role  in  the  characteristic 
of  flow  behavior.  It  will  be  seen  such  a  constraint  may  be  responsible  for  the  variation  or  redistribution 
of  the  density  and  for  the  special  vortex  pattern  developed  in  the  stratified  late  wake  behind  a  sphere. 
The  effect  on  the  density  distribution  by  the  presence  of  a  pressure  constraint  can  be  qualitatively  seen 


from  the  following  discussion. 
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Consider  a  line  vortex  with  its  axis  of  symmetry  coinciding  with  the  z-axis  of  a  cylindrical  coordi¬ 
nate  system  (r,  0,  z).  The  constraint  equation  for  the  present  flow  is 

(p„ril2)  +  ~  ip0g)  •  0  (I) 

dz  or 

where  p()  and  H  are  respectively  the  density  and  the  angular  velocity  of  the  flow.  Let  us  first  assume 
that  the  density  is  stratified  only  in  the  axial  direction,  i.e.,  p0  -  pa(:).  The  constraint  equation 
requires  the  angular  velocity  to  be  described  by 

n-a (2) 

where  a >(r)  is  an  arbitrary  function  of  the  radius.  The  steady-state  pressure  is  now  governed  by 

P(rtz)  -  J  ra)2(r)dr  -  g  J*  p0(z)dz .  (3) 

Now  we  consider  two  vortices,  one  on  top  of  the  other,  with  a  common  axis  and  a  common  boundary 

located  at  r  *  Z.  The  vortex  on  top  has  a  density  p\  and  an  angular  velocity  fl  \  while  the  one  at  the 

bottom  has  a  density  p2  and  an  angular  velocity  n  2-  The  pressures  in  each  individual  region  are, 

respectively, 

/»,-  /  rwfrfr  -  g  j  Pidz  (4) 

*  J  rml  dr  -  g  J  p2dz 

The  pressure  balance  condition  at  the  common  boundary  z  *  Z  requires  that 


This  implies 


or 


Pj(r.Z)  -  P2(r,Z) 

co f ( r )  *  cof(r) 


p](Z )  ft  f(r)  -  p2(Z)  Q  Hr)  -  constant  (5) 

In  other  words,  if  the  density  were  restricted  to  be  z-axis  dependent  only,  the  angular  velocity  would  be 

inversely  proportional  to  the  density  of  the  fluid.  Such  a  restriction  in  turn  implies  that,  for  a  statically 

stable  density  distribution  of  a  line  vortex,  the  rotational  velocity  should  be  large  on  top  and  small  at 

the  bottom.  This  is  apparently  not  the  vortex  pattern  developed  in  the  late  wake  behind  a  sphere.  The 

above  discussion  suggests  that,  in  order  to  generate  the  particular  vortex  pattern  behind  an  axisym- 
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metric  body,  the  steady-state  density  will  have  to  be  redistributed  from  its  original  axis-dependent  dis¬ 
tribution.  The  density  difference  between  the  wake  and  the  environment  would  be  even  sharper  as  a 
result  of  such  a  redistribution.  This  could  be  the  cause  for  the  existence  of  the  particular  vortex  struc¬ 
ture  far  down  stream  of  the  wake. 

To  understand  the  behavior  of  such  a  vortex  pattern  behind  an  axisymmetric  body,  one  needs  to 
investigate  the  stability  characteristics  of  a  general  class  of  flows  which  have  their  density  and  velocity 
components  varying  in  more  than  one  direction. 

A  classical  way  to  attack  problems  for  flow  stability  is  to  apply  the  well-known  normal  mode 
method  to  a  set  of  partial  differential  equations  governing  the  stability  characteristics  of  flows.  The 
method  is  very  effective  for  flow  profiles  with  one  independent  variable  since  periodic  solutions  are 
readily  admitted  to  the  partial  differential  equations  which  then  are  reduced  to  a  set  of  ordinary 
differential  equations.  For  flow  profiles  with  density  and  velocity  distributions  depending  on  more  than 
one  coordinate,  the  method  of  normal  modes  becomes  exhausted  because  the  partial  differential  equa¬ 
tions  governing  flow  stability  in  general  can  not  be  reduced  to  a  corresponding  set  of  ordinary 
differential  equations.  The  method  of  separation  of  variables  may  provide  a  way  to  solve  the  problem, 
however,  its  application  is  limited  to  only  a  few  special  cases. 

The  method  of  generalized  progressing  wave  expansion  can  circumvent  this  difficulty  and  provides 
us  with  a  tool  to  attack  stability  problems  for  flows  varying  along  more  than  one  coordinate.  In  a  recent 
paper  by  Eckhoff  and  Storesletten  (1978),  this  method  was  applied  to  a  class  of  inviscid  helical  gas  flows 
with  their  steady-state  profiles  depending  oniy  on  the  radius.  A  necessary  condition  for  stability  was 
derived  and  compared  with  the  existing  stability  criteria.  In  the  present  investigation,  a  general  class  of 
compressible  vortex  flows  with  their  steady  state  distributions  depending  on  the  radial  and  axial  coordi¬ 
nate  are  to  be  considered.  Dissipation  effects  due  to  viscosity  and  thermal  diffusivity  are  disregarded. 
Necessary  conditions  for  stability  of  the  flows  are  derived  using  the  method  of  generalized  progressing 
wave  expansion.  For  the  flows  to  be  stable,  it  is  necessary  that  they  be  stable  in  the  centrifugal  force 
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field  created  by  the  rotation  of  the  fluids,  in  the  gravitational  force  field  arisen  from  the  density  varia¬ 
tion  in  the  vertical  direction,  and  in  the  pressure  field  balanced  by  centrifugal  and  gravitational  forces. 
These  conditions  can  be  interpreted  by  a  kinetic  energy  approach  similar  to  the  one  used  by  Rayleigh 
(1916)  and  by  a  work  done  approach  based  on  the  movement  of  fluid  particles  in  the  centrifugal  and 
gravitational  force  fields. 


We  have  not  been  successful  in  obtaining  a  sufficient  condition  for  stability.  However,  it  is  sug¬ 
gested,  from  the  necessary  condition  to  be  derived,  that  a  sufficiency  condition  for  flows  with  the  steady 
distributions  varying  in  both  the  radial  and  axial  directions  should  satisfy  three  criteria:  one  in  the 
radial  direction  (Fung  and  Kurzweg,  1975;  Lalas,  1975),  one  in  the  vertical  direction  (the  classical 
Richardson  criterion  for  2D  parallel  flows),  and  a  third  one  constrained  by  the  pressure  relation  for  the 
variation  of  the  density  and  for  the  balance  of  the  two  force  fields. 

GOVERNING  EQUATIONS 


Consider  an  isentropic  vortex  motion  having  the  axis  of  symmetry  aligned  with  the  r-axis  of  a 
cylindrical  coordinate  (r,  9,  z)  and  with  the  direction  of  gravity.  The  fluid,  rotating  at  an  angular  velo¬ 
city  fl  (r,  z)  and  having  a  density  p„(r,  c),  is  assumed  to  be  compressible  but  inviscid.  To  satisfy  the 
pressure  balance  condition  anywhere  within  the  flow  regime,  the  constraint  equation  (1)  for  the  steady- 
state  must  be  satisfied.  Assume  the  perturbations  to  the  steady-state  flow  profile  are  small.  The  linear¬ 
ized  partial  differential  equations  governing  the  system  are 
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D  d  S 

where  —  —  —  +  a  —  and  c0  is  the  velocity  of  sound.  The  perturbations  of  velocity  in  r , 
Dt  ot  d9 

directions  are,  respectively,  u,  v  and  w  while  the  perturbations  of  the  density  and  pressure  are, 
tively,  p  and  p. 


Consider  the  following  transformation  for  the  density  and  pressure  perturbation  such  that 


p  -  —  (p  +  p) 

Co 

P  “  C0p0p. 

Equations  (6)  can  be  transformed  into  a  symmetric  hyperbolic  system  described  by 
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Following  the  method  of  generalized  progressing  wave  expansion  used  by  EckhofT  and  Storesletten 
(1978)  in  their  study  on  the  stability  of  helical  gas  flows,  we  obtain  the  characteristic  equation  for  the 
symmetric  hyperbolic  system  described  by  (10)  as  follows: 

det  i£{A|  +  f2A2  +  f  jA3  -  Xl|  -  0.  (9) 

The  characteristic  roots  of  the  determinant  are  given  by 


xj -  n*2 

X :  -  n  f  2  +  c<>V  (10) 

x 3 *  nf2-  c„v. 

Here  f  ,  are  components  for  the  orthonormai  eigenvectors  associated  with  the  eigenvalue  X  and 

v-  dr  +(i/r2  +i})U2  (11) 

The  eigenvalues  X2  and  X3  correspond  to  acoustic  solutions  and  will  not  be  discussed  any  further.  The 

eigenvalue  X  t  corresponds  to  gravity  waves  and  is  the  main  concern  in  this  paper.  The  corresponding 

ray  equations  for  gravity  waves  are 
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with  the  solutions  to  these  ray  equations  given  by 
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where  r„,  r)„,  2„  and  f  (Ac  —  1.2,3)  stand  for  the  initial  values  at  r  *  0.  The  amplitude  of  the  leading 
term  of  the  generalized  progressing  wave  expansion  for  the  gravity  waves  is  described  by 
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where  cr^.  are  scalar  functions  to  be  used  as  dependent  variables  in  the  transport  equations,  and  rK  are 


orthonormal  eigenvectors  defined  by 


r /  •  Ai  r, 


The  eigenvectors  are  selected  such  that 
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The  transport  equations  for  the  flow  under  consideration  are  given  by 
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Here  <r  =  {cr ,,  <r:,  <x:,}  and  A(/)  is  a  time  dependent  matrix  with  its  elements  described  by 
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0)2 
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The  transport  equations  are  valid  along  the  rays  described  by  Eqs.  (13).  The  stability  characteristics  of 
the  (low  governed  by  partial  differential  Eqs.  (6)  are  now  equivalent  to  those  of  ordinary  differential 
Eqs.  (18)  evaluated  along  the  rays.  The  flow  will  not  be  stable  unless  the  system  (18)  is  stable  for  all 
possible  choices  of  r0%  z0  and  £  k0. 

STABILITY  CHARACTERISTICS 


The  system  described  by  the  transport  Eq.  (18)  is,  in  general,  non-autonomous.  A  stable  asymp¬ 
totic  behavior  for  large  t  will  be  required  in  order  to  have  stability  for  the  system.  For  /—♦*>,  the 
non-vanishing  elements  of  matrix  A  are 


Frx 2  4-  [F:  4-  D*(rQ  )1  xy  4  —  (rfl  )y2 

dz 


a  23 : 


<*)\ 


D*(rQ.)  x2  - 


Fr~  (rQ) 

a: 


xy  -  F.y2 


•/A 

•/A 


L  iHi  x 2  _ 

2ft - — 

xy 

1 

c„ 

—  20  x2  +  — — —  .xv  —  y2 

C„  c„ 


/A 

■/A 


(21) 


9 


r 


Y  T  FUNG 

where  x  -  y  “  and  A  -  x2  +  y2.  The  characteristic  equation  for  the  eigenvalue  X  of  the 
or  or 

non-autonomous  system  described  by  Eq.  (21)  can  be  found  as 

X(X2  -h'l'/A)  -  0 


where 


* 


F,  — —  +  2ft  ft*(rft ) 

Cn 


X “  - 


Fr  x.  -  f.  -  2ft  4~  < m ) 

C„  ‘  C„  dz 


xy 


(22) 


-F.-*-y2.  (23) 


In  order  that  the  system  with  the  elements  described  by  Eqs.  (21)  shall  be  stable,  it  is  necessary  that  all 
the  characteristic  roots  in  Eq.  (22)  have  non-positive  real  parts,  i.e.. 


•*  4 

•  i 

A 

:<  ; 


♦  Ss  0. 

Considering  .vand  yas  two  independent  variables,  one  may  normalize  condition  (24)  into 
Fr  +  2ft  Z)*(rft )  —  F.  —  ± 

Cn  '  C„ 


(24) 


Fr  — -  +  2n  D*(rft )  —  F. — 

Cn  "  C„ 


Equation  (25)  will  be  satisfied  if 


+  4F. 


F, -^  +  2ft£>*(rft) 

Cn 


+  4 


Fr*- 

Co 


> 0 .  (25) 


Fr  +  2ft  Z)*(rft )  -  F.  —  >  0 

Cn  *  C„ 


and 


F.  -S- 

F,  +  2ft  ( rft) 

F,  X 

'  c„ 

c„ 

(26a) 


(26b) 


Since  the  right  hand  side  of  Eq.  (26b)  is  always  positive,  conditions  (26)  incorporated  with  Eqs.  (20) 
require 

12 


/•n2  ap0 

P»  dr 


+  — 


rft  ‘ 


>  0 


X  ^  + 

121 

X 

p0  a.* 

^0 

>  0 


(27a) 

(27b) 


and 


_g_  ap^  + 

X 

'1 

/•ft2  ap„ 

/•ft2 1' 

tl 

g  \ 

rft2  ap„  _ 

rft2  J‘ 

Pn  <)Z 

Cq 

P.»  Br  + 

rft2 

P,.  ®r 

Co 

(27c) 


10 


fi*i' '  ■  i  i'i *  tiajhAiii  mI  ■  ai 


NRL  MEMORANDUM  REPORT  4571 


Here  <t>  —  2 HD*(rfl)  reduces  to  the  classical  Rayleigh  discriminant  if  the  azimuthal  velocity  of  the 
flow  is  radius-dependent  only. 


For  uniform  rotation,  i.e.,  fl  —  constant,  no  criteria  can  be  observed  since  Eq.  (24)  will  become 
trivial.  Special  considerations  will  have  to  be  taken  in  order  to  obtain  a  general  stability  criterion  for 
every  possible  flow  profile  under  the  present  assumption.  For  dfl/dr  —  bFl/bz  -  0,  system  (18) 
becomes  autonomous,  and  the  ray  equations  (13)  for  autonomy  reduce  to 

k-  1,2,3  (28) 

The  eigenvalues  X  of  system  (18)  are  now  governed  by  the  equation 

AlX-’  +  tfAfft  +  ffo  +  f/o)!-  0  (29) 

where 


F,  — — -  +  4ft 2 

ii'o- 

If,  g  f.  rii: 

F:  + 

F,  —  —  F. 

Z  :o 

Co  '  Cit 

Co 

C,t  C0 

r  J 

In  order  to  have  system  (18)  be  autonomous  and  stable,  it  is  required  that  the  eigenvalues  have 
non-positive  real  parts,  i.e., 

*  >  0.  (30) 

Since  f  )0,  Zzo  and  £30  are  independent  components  of  the  orthonormal  eigenvectors,  inequality  (30)  for 

arbitrary  values  of  Z  io*  Z  :o  and  Z  so  is  equivalent  to 

F,  -  F.  >  0  Ola) 

"  c„ 

and 


Fr  +  4n 2 

£,Jo- 

F,  -  F. 

Co 

Cq  Cn 

Z  10  Z  30 


-  F;  if*  >  0 


(31b) 


Equation  (31b),  as  a  quadratic  equation  for  arbitrary  values  of  f  10  and  £30.  can  be  normalized  into 
the  following  two  inequalities 


F. 


+  4H 2 


-  F.  -Z-  >  0 

'  c0 


F  rn2 

r(l 2 

2 

F.  J- 

'  c„ 

g 

Co 
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(32a) 

(32b) 
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It  can  be  seen  immediately  Eqs.  Ola)  and  (32a)  will  be  satisfied  automatically  if  Eq.  (32b)  is  satisfied. 
To  have  stability  of  the  system  of  equations  (18)  for  autonomy,  it  is  therefore  required  that 


rft;  dp«  _ 

Pa  dr  c0 


-  J-  ^  +  Ufl  0. 

Po  OZ  CQ  j 


Equations  (27)  and  (33) ,  evaluated  along  the  rays  governed  by  Eqs.  (13),  represent  the  stability 
conditions  of  the  system  of  equations  (18)  for  both  non-autonomous  and  autonomous  cases,  respec¬ 
tively.  The  stability  characteristics  of  the  flow  governed  by  the  partial  differential  Eq  (6)  are  equivalent 
to  those  of  the  system  governed  by  the  ordinary  differential  Eq.  (18)  evaluated  along  the  rays.  It  can 
therefore  be  concluded  that  the  necessary  conditions  for  stability  of  the  compressible  stratified  flows 
under  consideration  require  the  inequalities 

N2  +  0  ^  0  (34a) 

N2  >  0  (34b) 

(Nr2+<t>)Nz2  3*  NrVF2  (34c) 

to  be  satisfied  anywhere  within  the  flow  regime.  Here  the  Brunt-Varsala  frequencies  are  defined  as 


r2  ra2  dp0  [  rQ 2 


Po  dr 


fsj}  _ _ S_ 


j_  dp_o_ 
Po  dz 


to  respectively  represent  the  density  variations  along  the  radial  and  axial  force  fields.  The  Froude 
number  is  defined  as 


to  represent  the  ratio  between  the  centrifugal  and  gravitational  forces.  Equation  (34a)  will  be  replaced 
by  S;  >  0  for  uniformly  rotating  flows.  It  is  worth  mentioning  that  Eqs.  (34)  actually  do  not  represent 
three  independent  conditions.  Condition  (34c)  implies  that  either  one  of  the  two  conditions  in  (34a) 
and  (34b)  must  be  automatically  satisfied  if  the  other  one  is  satisfied.  This  behavior  as  well  as  the  phy¬ 
sical  mechanism  behind  the  three  conditions  will  be  explained  in  the  following  section. 
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INTERPRETATION  OF  THE  RESULTS 

The  necessary  conditions  for  stability  previously  derived  can  be  explained  based  on  a  kinetic 
energy  approach  similar  to  the  one  used  by  Rayleigh  (1916)  and  on  a  work  done  approach  based  on  the 
movement  of  fluid  particles  in  the  centrifugal  and  gravitational  force  fields.  For  simplicity,  compressi¬ 
bility  effects,  in  general  destabilizing  the  flow  as  implied  by  conditions  (34),  will  be  ignored  in  the  fol¬ 
lowing  explanation  of  the  stability  conditions.  It  should  be  first  pointed  out  that,  for  the  inviscid  flow 
under  consideration,  the  principle  of  conservation  of  circulation  is  valid,  since  the  density  is  single¬ 
valued.  This  is  the  key  point  for  the  physical  arguments  of  the  stability  condition  to  be  immediately 
discussed. 


Consider  two  fluid  particles  originally  located  at  Q{  and  Qz  within  the  flow  regime  in  the  r-z  plane 
as  illustrated  in  Fig.  2.  The  particle  at  Q \  has  a  density  pa  and  a  velocity  v^  while  the  particle  at  Q:  has 

9  9 

a  density  p()  4*  5  p0  and  a  velocity  v0  4*  8v0.  Here  v0  *  r  ft  and  8  -  8r  —  +  8z  To  explain  the 

or  oz 

physical  meaning  of  the  necessary  conditions  for  stability,  we  first  use  the  energy  approach  by  consider¬ 
ing  the  variation  of  the  total  energy  as  a  result  of  a  perturbation  to  the  system.  In  the  steady-state,  the 
kinetic  and  potential  energy  of  the  two  particles  are  given  by 

K.E.  zr  (p o  V(>  +  (p„  +*P0)(va  +8va)2I 

2  (37) 

and 


P  E.  p0gz  4*  ip0  +SpQ)  g(z  +$r). 

When  the  two  particles  interchange  their  positions  as  a  result  of  a  perturbation,  the  kinetic  and  the 
potential  energy  of  the  perturbed  system  are 


K.E. 


and 


I 

jp<> 

rv„ 

2 

r  4-  8  r 

i  +  SPo) 


(r  4  8r)(v„  4*  8v„) 


P.E.  +8r)  4-  (p0  48p0)gz. 


(38) 
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Fig.  2  —  Coordinate  of  the  fluid  particles 


Here  the  conservation  of  circulation  has  been  applied  to  the  resultant  kinetic  energy  in  (38).  If  the  per¬ 
turbation  is  small,  the  stability  of  the  system  requires  that  the  leading  term  of  the  energy  variation  as  a 
result  of  interchange  of  the  two  particles  be  equal  or  greater  than  zero,  i.e.. 


rfl*  dp„ 
- r —  +  <t> 

(dr)1- 

Po  dr 

everywhere  within  the  flow  domain. 


dp0 

dr 


J_i_  , 

a  _  ^P  i: 

Po  OZ 


n:) 


(8r)(8r) 


~  (8r)2  >  0 
Po  dz 


(39) 


An  alternative  approach  one  may  use  to  observe  the  stability  characteristics  of  the  system  is  to 
examine  the  work  done  by  the  two  particles  in  the  centrifugal  and  gravitational  force  fields.  When  the 
two  particles  interchange  their  positions,  the  work  done  by  the  particle  originally  located  at  Qj  is 


(r  v„): 

( r  4 -  Sr)’ 


+ PoS&: 


(40) 
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while  the  work  done  by  the  particle  originally  located  at  Q:  is 
Po  +Sp0 


w,  > 


(v„  +  8vu)2  [(r  +  8r)(v„  +  8v,,))2 


(r  +  8r) 


Sr  -  <P„  +Sp„)#8z. 


(41) 


Condition  (39)  can  be  reached  following  the  argument  that  the  stability  of  the  system  requires  the  lead¬ 
ing  terms  of  the  total  work  done  by  the  interchange  of  the  two  particles  to  be  non-negative,  i.e.. 


W, 


+  Wz  -  (p„r2v2]8r  +8<p0g)8:  >0 

*  J 


(42) 


The  physical  meaning  of  the  necessary  condition  for  stability  derived  by  the  method  of  general¬ 
ized  progressing  wave  expansion  can  be  seen  immediately  if  one  compares  Eq.  (24)  with  Eq.  (39)  or 
(42).  Equation  (24)  for  incompressible  fluids  reduces  to 

ll  an  Kan  1  f  ~  An  II  an  1^ 

>0.  (43) 


rft  ’ 

fr-  ,A 

an 

2 

g  1 

f>0 

dr  +*\ 

dr 

Po  dr  p„ 

an 

an 

g  &P° 

an 

1  dr 

dz 

Po  dz 

dz 

Equations  (39)  and  (43)  are  identical  if  one  considers  5rand5zas  independent  variables 

or  dz 

and  match  with  and  with  Sz.  If  the  variation  in  the  axial  direction  is  restrained ,  Eqs.  (39) 
or  oz 

and  (43)  reduce  to  the  well-known  Rayleigh-Synge  criterion  (Synge  1933)  which  is  a  requirement  for 
centrifugal  stability.  Also  it  can  be  easily  shown,  by  the  same  integral  method  used  by  Fung  (1980), 
that  the  first  term  in  Eqs.  (39)  or  (43)  is  in  fact  a  differential  representation  of  a  stable  centrifugal  force 
field.  This  mechanism  is  reflected  in  Eq.  (34a)  saying  that  the  steady  state  of  the  flow  should  be  stable 
in  the  radial  direction.  Two  parts  are  involved  in  this  first  term.  The  first  part  is  the  variation  of 
density  in  the  centrifugal  force  field.  The  second  part  is  the  Rayleigh  discriminant  representing  a  varia¬ 
tion  of  the  centrifugal  balance  as  a  result  constrained  by  the  conservation  of  circulation. 


If  the  variation  in  the  radial  direction  is  suppressed,  Eqs.  (39)  and  (43)  reduce  to  a  condition 
representing  the  variation  of  density  in  the  gravitational  force  field.  This  mechanism  is  reflected  in 
(34b)  saying  that  the  steady  state  of  density  should  be  statically  stable  in  the  axial  direction. 

The  second  term  in  Eqs.  (39)  and  (43)  represents  the  cross  correlation  between  the  radial  and 
axial  variations  of  the  density  in  the  centrifugal  and  gravitational  force  fields.  This  variation  is  reflected 


15 


Y.  T.  FUNG 


in  condition  (34c)  with  a  mechanism  to  be  immediately  revealed.  For  incompressible  flows,  Eq.  (34c) 
reduces  to 


4> 

-  i  +  — 

rQ2  df>o 
P 0  ^ 

Equation  (44)  represents  a  requirement  for  stability  imposed  on  the  density  variations  in  the  radial  and 
axial  directions.  The  second  term  inside  the  bracket  in  the  above  equation  is  the  Rayleigh  discriminant 
normalized  by  the  centrifugal  force  as  a  result  of  the  density  variation  in  the  radial  direction.  It  is  a 
ratio  of  the  variation  of  centrifugal  forces  between  that  resulting  from  the  conservation  of  circulation 
and  that  arising  from  the  variation  of  density  in  the  centrifugal  force  field.  For  fixed  values  of  the 
Froude  number,  the  density  variation  in  the  radial  direction,  as  compared  to  that  in  the  axial  direction, 
plays  an  opposite  role  in  flow  stability.  While  large  density  gradients  (positive)  in  the  axial  direction 
stabilize  the  flow  as  in  the  case  of  two  dimensional  stratified  fluids  in  a  gravitational  force  field,  large 
density  gradients  (positive)  in  the  radial  direction  destabilize  the  flow.  This  is  opposite  to  the  role 
played  by  density  gradients  in  radius-dependent  rotating  flows,  and  seems  to  be  implausible  at  the  first 
look.  However,  this  is  also  an  interesting  point  that  will  reveal  the  physical  mechanism  carried  by  Eqs. 
(34c)  and  (44).  For  potential  flows,  both  Eqs.  (1)  and  (44)  reduce  to 


ap.  f 

dr 


(44) 


dpo  /dp,, 

dr/  dz 


ril2 

g 


(45) 


saying  that  the  ratio  between  the  density  gradient  in  the  radial  direction  and  that  in  the  axial  direction 


should  be  compatible  with  the  ratio  between  the  centrifugal  force  field  and  the  gravitational  force  field. 
And  it  is  also  the  constraint  condition  for  the  pressure  balance  requires  that  large  density  gradients 
(positive)  in  the  radial  direction  result  in  small  density  gradients  in  the  axial  direction,  destabilizing  the 
flow. 


Based  on  the  arguments  just  presented  for  the  physical  mechanisms  of  the  necessary  conditions, 
we  conclude  that  Eqs.  (34)  represent  a  generalized  state  of  statically  stable  profiles  for  the  steady  flow. 
To  secure  stability  for  the  basic  flow,  it  is  necessary  that  the  steady-state  distribution  satisfy  the  radial 
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force  balance  condition,  the  axial  force  balance  condition,  and  a  pressure  balance  condition  constraining 
the  variations  of  density  in  both  the  centrifugal  and  gravitational  force  fields.  As  a  result  of  the  third 
constraint,  Eqs.  (34)  do  not  represent  three  independent  conditions.  Either  one  of  the  two  conditions 
in  (34a)  and  (34b)  will  have  to  be  automatically  satisfied  if  the  other  condition  and  Eq.  (34c)  are 
fulfilled. 

CONCLUSIONS  AND  DISCUSSIONS 

Necessary  conditions  of  stability  for  an  isentropic  compressible  vortex  flow,  with  the  density  being 
stratified  in  both  the  radial  and  axial  directions,  have  been  derived  using  the  method  of  generalized  pro¬ 
gressing  wave  expansion.  This  method  transforms  a  set  of  partial  differential  equations  into  a  set  of 
ordinary  differential  equations  with  equivalent  stability  characteristics  evaluated  along  the  rays.  Thus  it 
provides  one  with  a  powerful  tool  to  attack  a  class  of  steady  flows  with  more  than  one  independent  vari¬ 
able,  in  which  the  method  of  normal  modes  is  exhausted.  The  necessary  conditions  derived  here 
represent  a  generalized  state  of  statically  stable  distribution  for  the  steady  flow.  They  require  the  flow 
to  satisfy  the  centrifugal  force  balance  condition,  the  gravitational  force  balance  condition,  and  the  pres¬ 
sure  balance  condition  restraining  the  variation  of  densities  and  forces  in  both  force  fields. 

It  is  shown  in  the  Appendix,  by  using  the  Buckingham  n  theorem,  that  the  general  class  of  vortex 
flows  being  considered  possesses  three  non-dimensional  numbers.  They  are  the  Mach  number,  the 
Froude  number,  and  the  Richardson  number  (the  Reynolds  number  will  appear  if  viscous  effects  are 
considered).  The  first  two  are  brought  out  by  the  necessary  conditions.  The  last  one  is  obviously 
embedded  in  a  sufficiency  condition  for  stability  yet  to  be  derived.  One  way  to  approach  such  a 
sufficiency  condition  is  to  construct  a  Liapunov  function  for  the  system  of  equation  (14)  using 
Liapunov’s  direct  method.  We  have  not  been  successful  in  such  an  approach. 

Even  though  it  has  not  been  proved,  several  implications  concerning  the  sufficiency  condition  or  a 
Richardson  criterion  can  be  observed  from  the  necessary  conditions.  In  view  of  Eqs.  (34),  the 
sufficiency  condition  should  also  consist  of  three  criteria:  one  in  the  radial  direction  (Fung  &  Kurzwcg, 
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1975;  Lalas,  1975),  one  in  the  vertical  direction  (the  classical  Richardson  criterion  for  2D  parallel 
flows),  and  a  third  one  imposed  by  the  pressure  constraint  equation.  Such  a  conjecture  is  based  on  the 
necessary  and  sufficient  conditions  for  stability  so  far  derived  for  two-dimensional  stratified  parallel 
flows  and  radially  stratified  rotating  flows.  For  two-dimensional  incompressible  flows  with  densities 
stratified  in  z-direction,  the  necessary  condition  for  stability,  or  ^statically  stable"  as  usually  called,  is 


the  corresponding  sufficiency  condition  is 


(46) 


dpQ 

dz 


dU 


dz 


(47) 


where  U  is  the  parallel  velocity  perpendicular  to  the  direction  of  gravity.  For  incompressible  rotating 
flows  with  densities  and  velocities  varying  only  in  the  radial  direction,  the  necessary  condition  for  stabil¬ 
ity  is 


or 


4> 


r ft2  dPo 
Pa  dr 


+  4> 


>  0 


(48a) 


dp  „ 
dr 


>  0 


(48b) 


for  uniformly  rotating  flows.  Here  4>  is  the  well-known  Rayleigh-Synge  discriminant  and  Eq.  (48a) 
becomes  the  necessary  and  sufficient  condition  for  stability  if  only  axisymmetric  disturbances  are  con¬ 
sidered.  The  corresponding  sufficient  conditions  for  stability  of  the  flow  subject  to  azimuthal  distur¬ 
bances  is 


<D 


_  >  i 

( rDCl  +  rCl)2  '  4 


(49a) 


or  in  an  alternative  form 


rfl  * 


dp,> 

dr 


>  i 


(49b) 


pJrDtl)2  4 

Based  on  conditions  (30)  and  on  the  stability  criteria  just  discussed,  we  may  propose  the  following 
sufficient  conditions  for  stability  of  the  genera)  class  of  flows  being  considered. 


18 


SRL  MEMORANDUM  REPORT  4571 


If  the  Brunt-Vaisala  frequencies  are  defined  as  in  Eqs.  (35)  to  respectively  represent  the  density 
variations  in  the  radial  and  axial  directions,  the  appropriate  local  Richardson  numbers  in  the 
corresponding  direction  may  be  defined  as: 


and 


1 


Nr2+<t> 


r^P~  +  4fl 

or 


(50) 


J: 


N? 


(51) 


The  plausible  sufficiency  conditions  for  a  system  subject  to  azimuthal  disturbances,  if  proved,  would 
have  the  forms 


J  >  1 

J'  *  4 

(52a) 

J.  >  \ 

(52b) 

4 

fUr.J:)  >  0 

(52c) 

where  f{JrJ:)  is  an  arbitrary  function  of  the  two  Richardson  numbers.  Equation  (52a)  is  a  Richardson 
criterion  in  the  radial  direction  and  can  also  be  written  as 


N2 


dn_ 

dr 


(53) 


Equation  (52b)  is  a  Richardson  criterion  in  the  axial  direction  while  Eq.  (52c)  is  a  pressure  balance 
requirement  constraining  the  two  Richardson  criteria  in  both  directions.  As  implied  by  the  necessary 
conditions  for  stability  in  Eqs.  (34),  the  criteria  described  by  Eqs.  (52)  would  not  represent  three 
independent  conditions.  Either  one  of  the  two  conditions  in  (52a)  or  (52b)  would  automatically  be 
satisfied  if  the  other  condition  and  Eq.  (52c)  are  fulfilled.  Furthermore,  the  necessary  conditions  will 
serve  as  a  prerequisite  for  the  derivation  of  the  sufficiency  conditions  yet  to  be  identified. 


An  understanding  of  the  stability  criteria  for  the  flow  under  consideration  is  essential  to  the  vor¬ 
tex  motion  in  the  late  wake  behind  a  towed  axisymmetric  body.  The  existence  of  the  vertically  oriented 
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vortex  structure  far  downstream  of  the  wake  implies  that  the  vortex  motion  satisfies  stability  criteria 
such  as  those  given  by  Eqs.  (52).  The  stability  criteria,  if  derived,  will  provide  us  with  an  insight  into 
the  flow  characteristics  that  may  be  used  to  predict  the  evolution  and  breakdown  of  the  vertically 
oriented  vortex  motion. 
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APPENDIX 


The  dimensionless  numbers  for  inviscid  and  compressible  vortex  motion  in  the  stratified  wake 
behind  an  axisymmetric  body  are  derived  as  follows.  The  independent  dimensions  for  the  flows  are: 

M  =  mass 

L  **  characteristic  length 
T  *  time. 

The  independent  variables  are: 

p()  =  density  (M/L3) 

V  —  velocity  (L/T) 

D  *  characteristic  length  of  the  body  (L) 
g  =*  gravitational  constant  (M/TL) 
cr,  =  velocity  of  sound  (L/T) 

=  Brunt-Vaisala  frequency  (1/T2) 

The  Buckingham  rr  theorem  states  that  the  number  of  independent  dimensionless  numbers  is  equal  to 
the  number  of  independent  variables  less  the  number  of  independent  dimensions,  i.e.,  6-3*3. 
These  three  independent  dimensionless  numbers  are  constructed  in  the  following  way: 


Dimensionless  number 
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a  :  —  —  l 
a  t  *  0 
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and 
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Dimensionless  number 
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